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ABSTRACT 

Four  right-hand  orthogonal  coordinate  systems  are  defined  and  the  trans¬ 
formations  between  these  systems  are  developed.  They  are:  a  local 
topocentric  launcher  system,  an  earth-centered  inertial  system,  a  body 
axis  systeip,  and  an  aeroballistic  axis  system. 

The  determination  of  time  derivatives  of  direction  cosines  is  used  to 
transform  information  between  the  inertial  and  body  systems.  Euler 
angles  are  used  to  transform  data  between  the  inertial  and  aeroballis¬ 
tic  axis  systems.  Strengths  and  weaknesses  of  these  transforms  are 
then  compared.  Finally,  a  comparison  of  results  of  trajectory  simu¬ 
lations  is  presented. 
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INTRODUCTION 

There  has  been  and  still  is  considerable  interest  in  the  simulation  of 
unguided  rocket  trajectories.  Several  mathematical  models  [l,  2D  have 
been  developed  for  this  purpose.  The  Atmospheric  Sciences  Laboratory 
has  developed  a  six-degree-of- freedom  ballistic  simulation  model  [3D. 
This  model  uses  three  right-hand  orthogonal  coordinate  systems  -  the 
launcher,  the  inertial,  and  the  "body"  coordinate  systems.  The  output 
from  the  simulations  is  referenced  to  the  launcher  coordinate  system. 
Aerodynamic  forces  and  moments  are  computed  in  the  body  system,  and 
the  inertial  system  is  used  for  the  development  and  integration  of  the 
equations  of  motion.  Information  is  trens formed  from  the  inertial  to 
body  system  via  direction  cosines.  Due  to  computational  inaccuracies, 
however,  it  is  possible  here  to  obtain  a  nonorthonormal  transformation. 
Hence  it  was  decided  to  determine  what  would  happen  if  fhe  body  axis 
system  were  replaced  by  an  aerobal I istic  axis  system  and  if  Euler  angles 
were  used  in  the  transformation  of  information  from  the  inertial  to  the 
aerobal I istic  axis  system.  This  report  compares  the  limitations  and 
capabilities  of  these  two  systems  and  their  results  in  theoretical  tra¬ 
jectory  simulations. 


INERTIAL  AND  LAUNCHER  COORDINATE  SYSTEMS 

Two  right-hand  orthogonal  coordinate  systems  will  be  defined.  (See  Fig¬ 
ure  I.)  The  launcher  coordinate  system  (denoted  X',  Y',  Z')  has  its 
origin  at  the  launcher  and  rotates  with  the  earth.  The  positive  X'-axis 
points  east,  the  positive  Y'=axis  points  north,  and  the  positive  Z'-axis 
points  outward  along  the  radius  vector  from  the  center  of  the  earth. 

The  inertial  coordinate  system  (X,  Y,  Z)  has  its  origin  at  the  center  of 
the  earth  and  is  oriented  so  that  the  X-Y  plane  lies  in  the  earth's 
equatorial  plane  with  the  positive  Y-axis  pointing  initially  through 
the  longitude  of  the  launcher.  The  Z-axis  is  coincident  with  the  earth's 
axis  of  rotation  and  positive  toward  the  North  Pole.  This  system  does 
not  rotate  with  the  earth. 

The  mathematical  simulation  model  considers  the  earth  to  be  an  oblate 
spheroid  with  an  equatorial  radius  of  6,378,423  meters  and  an  eccentricity 
of  .00672267.  The  transformation  from  the  launcher  system  to  the  iner¬ 
tial  system  is,  in  matrix  form. 
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where  X  is  the  geodetic  latitude  of  the  launcher,  u>  is  the  earth's  ro¬ 
tational  speed,  Rq  is  the  distance  from  the  center  of  the  earth  to  the 
launcher,  and  t  is  time  after  launch.  The  preceding  transformation 
and  its  inverse  are  derived  in  [23. 


THE  BODY  COORDINATE  SYSTEM 

The  "body"  coordinate  system  (x,y,z)  has  its  origin  at  the  center  of 
gravity  of  the  rocket.  The  x-axis  coincides  with  the  longitudinal  axis 
of  the  rocket  and  is  positive  toward  the  nose.  The  positions  of  the  y. 
and  z-axes  are  determined  by  the  rocket's  motion.  The  initial  positions 
of  these  axes  are  defined  as  follows:  Let  0  be  the  angle  between  the 
x-axis  and  the  positive  Z'-axis  measured  from  the  Z'-axis.  The  y-axis 
lies  in  the  X-Y  plane  and  is  positive  in  the  direction  of  positive  0. 

(See  Figure  I . ) 

The  transformation  from  the  inertial  to  the  body  coordinate  system  via 
direction  cosines  is  obtained  as  follows:  Let  (£|,  I2,  £3),  (m|,  m2,  m3), 
and  ( n | ,  n2,  n3>  be  the  respective  direction  cosines  of  the  body  system 
with  respect  to  the  inertial  system  [43.  Then 
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where  Kj  denotes  a  unit  vector  along  the  i  axis.  Since  this  is  an 
orthogonal  system,  its  inverse  is  its  transpose.  Thus, 

Kv  =  £.  K  +  m.  K  t  n.  K  .  (2) 

X  lx  I  y  I  z 

Differentiating  (2)  with  respect  to  time  yields 

0  =  £.  K  +  l.  K  +  m.  K  +  m.K  +  n.  K  +  n.K  .  (3) 

lx  lx  ly  ly  Iz  Iz 


Let  w  define  the  angular  velocity  of  the  body  system  relative  to  the 
inertial  system.  Then 

w  =  pK  +  qK  +  r« 

K  x  M  y  z 

where  n,  q,  r  are  the  x,  y,  z  components  of  w.  Hence 

K  =  K  r  -  K  q,  (4) 

x  y  z  ^ 


3 


and 


(4) 


r. 


K  =  K  q  -  K  p. 
z  x^  Y 

Substituting  (4)  into  (3)  yields 

0  =  K  (t.  -  rm.  +  qn . )  +  K  (m.  -  pn .  +  rl, ) 
x  I  I  ^  I  y  I  K  I  I 

+  K^(n |  -  q£|  +  pm | ). . 

Thus, 

£|  =  rmj-qnj 

m  j  =  pn|-r£|  (5) 

n (  =  q^j-pmj  , 


Similarly,  one  obtains 


t  = 


r  m.  -  qn. 


mj  =  p  n.  -  r£j 

nj  =  ql.  -  p  m. 


i  =  I,  2,  3. 


(6) 


Integration  of  the  nine  differential  equations  in  (6)  yields  the  desired 
direction  cosine  transformation. 


It  is  easy  to  see  that  the  aforementioned  numerical  integration  can  give 
rise  to  a  nonorthonormal  system  due  to  computational  inaccuracies.  Hence 
the  information  transformed  by  this  system  will  contain  errors  which 
will  accumulate  throughout  the  simulation  and  can  cause  inaccurate  tra¬ 
jectory  determinations.  Hence  it  is  necessary  to  somehow  force  this 
system  to  be  orthonormal.  This  can  be  accomplished  as  follows. 


Fi rst  set 


l  j 

a2  +  i2  +  i2) 1/2 


2.  3. 


Then  the  dot  product  of 

i  =  (1,,  i2,  i3) 
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and 

m  “  ( m j ,  m-,,  m^) 
is  taken,  call  It  DOT.  Then  set 

m(  =  m,  -  £j*D0T,  I  =  I,  2,  3. 


Final ly  redefining 

m j  =  m./(m(  +  rr^  +  )'  ,  i  =  I,  2,  3 

and  taking  the  vector  cross  product  of  £  and  m,  yields 
n  ( n  | ,  1^2)  n^) . 

This  system  is  now  orthonormal  and  is  the  system  used  to  transform  in¬ 
formation  between  the  inertial  and  body  axis  system;. 

Several  things  should  be  noted.  Only  the  differential  equations  for  the 
£  and  m  vectors  need  be  used.  Thus  there  are  but  six  differential  equa¬ 
tions  in  (6)  to  be  integrated.  Further,  the  m  vector  is  modified  to 
force  it  to  be  orthogonal  to  £  and  n  is  completely  determined  by  £  and  m. 
It  should  be  noted  that  it  would  have  been  just  as  easy  to  have  started 
with  the  m  or  n  vector  instead  of  £  and  thus  to  obtain  an  analogous,  but 
somewhat  different  system.  The  above  procedure  contains  artificiality 
which  leads  one  to  suspect  the  result  obtained  through  it.  To  alleviate 
this  problem,  it  was  decided  to  use  Euler  angles  and  an  aerobal I istic 
axis  system. 


EULER  ANGLES 

Let  0  be  the  angle  of  rotation  about  the  Z-axis  from  the  X  to  the  Y-axis. 
Then  from  Figure  2  one  sees  that 


Let  be  the  angle  of  rotation  about  the  yj-axis  from  the  x  to  the  Z-axis. 
It  follows  from  Figure  3  that 
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Rotation  About  Z-Axis 

i'igure  2 
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(8) 


Let  ip  be  the  angle  of  rotation  about  the  X2~axis  from  y ^  to  Z2-  From 
Figure  4  it  follows  that' 
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Combining  (7),  (8),  and  (9),  one  obtains 


cosij)  cos£  cos$  sin£  sinY  -sin  <p  cosy 
sin<f>  cos£  sin$  sin£  sinY  +cos  <f>  cosY 
-sine  cos  e  sinY 


cos$  sine  cosY  +simj>  sinY 
sint  sine  cosY  -cost  sinY 
cos  e  cos  Y 


( 10) 


This  is  the  transformation  from  the  body  to  the  inertial  system  using 
Euler  angles.  Equation  (10)  contains  many  trigonometric  functions  and 
thus  takes  considerable  time  to  evaluate  on  a  computer.  Unlike  the 
direction  cosine  matrix,  no  artificiality  is  used. 

Let  w  be  defined  as  before.  Then,  the  body  system  Euler  rate  matrices 
become 
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This  can  be  seen  by  examining  Figure  5. 


(II) 


AEROBALLISTIC  AXIS  SYSTEM 


An  aerobal I istic  axis  system  is  one  such  that  the  y  and  z  axes  are  not 
allowed  to  rotate  about  x.  Thus  is  a  constant,  ipQ,  ip  =  0,  and 
p  =  / p  dt  where  p  is  redefined  as  the  physical  rolling  rate  of  the  rocket. 
The  t ran s formal  ion  from  the  body  to  the  inertial  system  by  Euler  angles 
is  still  valid  although  ip  =  ipQ.  The  aerobal  I  istic  system  Euler  rate 
matrices  become 
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From  equations  (II)  or  (12)  it  is  easily  seen  that  the  Euler  rate  matrices 
can  become  singular  at  C  =  ir/2  or  its  inverse  matrix,  A(A),  to  be  unde¬ 
fined  at  ?  =  "/2.  Hence  one  must  be  extremely  careful  in  the  use  of 
such  a  system.  Also,  note  that  the  choice  of  the  aerobal I istic,  i.e., 
nonrolling,  coordinate  system  reduces  the  number  of  equations  to  be  in¬ 
tegrated  and  the  number  of  trigonometric  functions  to  be  evaluated,  since 
ip  =  0  for  all  t.  The  previously  discussed  definitions  for  the  Euler 
angles  were  chosen  to  avoid  singularities  of  the  Euler  rate  matrices 
for  White  Sands  Missile  Range  usage. 


ANALYSIS 

The  equations  of  motion  of  an  ungulued  rocket  as  developed  in  C2J  were 
used  in  each  of  the  two  following  nstances.  First  the  inertial, 
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launcher  and  body  coordinate  systems  were  used  in  one  of  the  mathema¬ 
tical  models  in  conjunction  with  the  direction  cosine  transformation. 

In  the  second  case  the  inertial,  launcher,  and  aerobal I i stic  axis  sys¬ 
tems  were  used  with  Euler  angles.  These  two  cases  will  be  referred  to 
as  cases  I  and  2  in  that  which  follows. 

Several  trajectories  were  simulated  for  both  cases,  and  their  results 
in  computational  time,  burn-out  positions  and  velocities,  and  impacts 
were  compared.  The  only  differences  between  these  simulations  were  the 
coordinate  systems  used  and  the  transformations  between  them. 

The  Aerobee  350  rocket  was  used  for  this  comparison.  This  rocket  is  a 
two-stage  unguided  high-altitude  probe  which  is  fired  near  the  vertical. 
Since  this  rocket  burns  out  and  then  coasts  out  of  the  atmosphere,  the 
equations  of  motion  of  a  rocket  were  reduced  to  three  degrees  of  freedom 
after  the  rocket  leaves  the  atmosphere  (61  km).  This  procedure  results 
in  no  loss  of  accuracy  and  was  followed  in  both  cases.  A  constant 
launcher  setting  was  assumed  throughout  which  resulted  in  a  peak  altitude 
of  ^487  km  and  impact  range  of  ^204  km. 

Rocket  trajectories  were  simulated  for  both  cases.  These  simulations 
were  made  with  and  without  wind  and  with  and  without  the  rocket  rolling. 
All  computations  were  performed  on  the  Univac  1108  computer. 

The  burn-out  positions,  burn-out  velocities,  and  impact  points  varied 
only  slightly  between  the  two  cases.  There  was,  however,  a  large  dif¬ 
ference  in  their  computational  time  on  the  computer  (Table  I). 


TABLE  I 


COMPUTATION  TIME 

(SEC)  FOR 

AEROBEE  350 

1  .D. 

CASE  1 

CASE  2 

No  wind.  No  roll 

45.7 

140.8 

No  w i nd,  Rol 1 

55.3 

155.7 

10  MPH  East,  No  Roll 

51.7 

222.8 

20  MPH  East,  Rol 1 

55.4 

248.9 

These  discrepancies  in  computer  running  time  can  be  explained  by  the 
large  number  of  evaluations  of  trigonometric  functions  needed  for  the 
desired  transformation  to  and  from  the  aeroballistic  axis  system  in 
Case  2.  On  the  average,  Case  2  required  four  times  as  much  computation¬ 
al  time  per  trajectory  simulation  as  Case  I. 
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It  should  be  noted,  however,  that  this  was  accomplished  on  a  slowly 
rolling  unguided  rocket.  For  such  rockets.  Case  I  is  vastly  superior 
to  Case  2  with  respect  to  computational  time  and  yet  seems  to  maintain 
acceptable  accuracy. 

Suppose  now,  one  wishes  to  simulate  the  trajectory  of  an  artillery  shell 
or  a  bullet  (something  with  a  rapid  roll  rate).  A  natural  question 
arises.  Are  Case  !  type  transformations  still  superior  to  those  of 
Case  2?  The  author  feels  that,  due  to  artificiality  in  the  determina¬ 
tion  of  the  direction  cosines,  Case  I  methods  can  lead  to  problems  and 
that  the  use  of  Euler  angles  will  prove  superior.  This  question  is  cur¬ 
rently  being  investigated.  Another  method  of  transforming  data  between 
two  coordinate  systems  is  by  the  use  of  quaternions.  This  method  also 
leads  to  an  orthonormality  problem  and  thus  was  not  considered. 


CONCLUSIONS 

In  comparing  the  two  cases,  several  strong  and  weak  points  naturally 
appear.  First,  in  Case  I,  due  to  the  necessity  of  the  numerical  in¬ 
tegration  of  the  derivatives  of  the  direction  cosines,  a  nonorthonormal 
system  can  arise.  The  artificial  forcing  of  the  direction  cosines  to 
be  orthonormal  leads  one  to  suspect  any  results  obtained  therefrom. 
However,  for  a  slowly  rolling  rocket,  the  method  of  Case  I  results  in 
a  quite  accurate  determination  (when  compared  to  other  methods)  of 
the  theoretical  rocket  flight  path  and  a  fast  simulation  time.  For 
rapidly  spinning  shells,  this  induced  artificiality  can  be  a  potential 
problem. 

With  regard  to  aercbal I istic  axis  and  Euler  angles,  one  must  contend 
with  the  evaluations  of  several  trigonometric  functions  to  obtain  the 
described  transformation.  This  case  has  none  of  the  artificiality  of 
the  first  case.  Near  the  singularity  of  the  Euler  rate  matrix,  many 
problems  can  and  do  arise.  Also,  the  method  is  slower  than  that  of 
direction  cosines  and  hence  requires  more  computational  time. 
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